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Процесс вытеснения нефти из пласта водой имеет ряд особенностей, 
обусловленных неоднородностью строения пласта и различием свойств 
вытесняющей и вытесняемой жидкостей. 

Простейшей моделью неоднородного пласта является пласт, 
состоящий из изолированных друг от друга прослоев, отличающихся 
различной проницаемостью. | 
р 7 Характерная особенность этой модели состоит в том, что каждый 


прослой является однородным, поэтому вытеснение нефти водой в каждом 


прослое подчиняется закономерностям для однородного пласта (рис. 1) [1]. 


» 


Схема движения фронта воды в слоистом пласте 


СД, ; 
| - нефть; 2 - смесь нефти с водой. 


Рис. | 


Неоднородность пласта по мощности (слоистость и расчлененность) 
обусловливает послойное вытеснение нефти и обводнение пласта, поэтому 
высокопроницаемые слои и пропластки оказываются обводненными 
раньше, чем менее проницаемые слои (рис. 2). Зональная неоднородность 
и его линзовидность обусловливают неравномерное обводнение пласта по“ 
площади. Наиболее проницаемые зоны и линзы пласта обводняются 
в первую очередь, на слабопроницаемых участках наблюдается очень 


медленное продвижение воды. 
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Типы перемещения водонефтяного контакта 


6 пластовых залежах нефти 


СЭ чальный ВНК 


8 


стена ___-- 2 


Пербонача ти ВНК 


1 
А“ Е 
а - нефть; 6 - вода, 1,2,3 - положения водонефтяного контакта. 


Рис. 2 


Картина движения воды в неоднородном пласте еще более 
усложняется, когда пласт скважинами вскрывается не полностью, 
а избирательно. 

В этом случае вода продвигается преимущественно по пропласткам, 
находящимся под воздействием закачиваемой воды и вскрытым 


эксплуатационными скважинами. 


1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ 
ДВИЖЕНИЯ ФРОНТА ВОДЫ В ПЛАСТЕ 


Движение водонефтяного контакта обусловлено перепадом давлений 
между нагнетательной и добывающей скважинами. Распределение 


давления по вертикальному срезу пласта может быть найдено из 


уравнения [3]: 
В.В = аіу(олр)- №, (1) 


где Дд. - коэффициент сжимаемости коллектора; 
р - давление; 


В - ширина коллектора; 


с - функция гидропроводности: 


ба ЕВ ИВ, | (2) 
ту е Ин Ив 
где К,К„„Ка - абсолютная проницаемость и относительные фазовые 
проницаемости нефти и воды соответственно; 
Др - градиент давления; 
М - функция, моделирующая работу скважин. 

Относительные фазовые проницаемости представляют собой 
однознанные — безразмерные функции насыщенности, Значения 
относительных фазовых проницаемостей в зависимости от насыщенностей 
определяются экспериментально и обычно аппроксимируются степенными 
зависимостями, причем наилучшее совпадение фактических и расчетных 
данных достигается при использовании квадратичных зависимостей: 

Кв = (58 Св) (1–а-(1–5р)) 
Е М Кн = (1-58 тен) (1-6-5) 


Ка 


(3) 


от къ 


где Св - насыщенность связанной ВОДЫ, 


Сн - остаточная нефтенасыщенность; 


- 


а, ђ - аппроксимационные коэффициенты, зависящие от свойств - 


фильтрующихся жидкостей и скелета коллектора. 


Движение жидкостей (нефть и вода) подчиняется закону Дарси: 


о 
У = – Ар, 4 
4 (4) 


где У - вектор скорости, тогда изменение насыщенности компонента 
может быть определено следующим образом: 


| да. 1 (е 
ЕЕЕ 5 


а нефтенасыщенность равна 


~ и Фа. Зи <1-5в. (6) 
2. ЧИСЛЕННАЯ РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ 


Полученная система уравнений может быть решена только численным 


методом. Уравнения в частных производных заменяются их конечно- 


разностными аналогами. 

| Для случая прямоугольной декартовой системы координат 

исследуемая область разбивается пох и у сеткой ячеек (плоский случай): 
х.12=(@+12)Ах, Ах>0, #=0,1,..,№; (М +1)Ах =. 


с Ура: = (7+1 2), Аг>0, ј=04,...,М; (М += а 


где Ах, Ах - размеры шагов сетки, 
№ М — число ячеек в направлении х, 2 соответственно (точка 
с координатами (Г, /) совпадает с центром ячейки). : 


Здесь используется «шахматная» сетка, т.е. координаты сеточнъх 


функций разнесены в пространстве. Это дает возможность наглядно’ 


интерпретировать каждую ячейку как элемент объема среды, который 


характеризуется рассчитываемыми в его центре давлением Рур» а также 


насъщенностями компонентов. Знание же нормальной составляющей 


пож ке: 


ПОСТАЛА А 
С Ди 


е 


+ 


вектора скорости на границе ячейки позволяет непосредственно вычислять 
поток вещества через эту границу. 

В качестве конечно-разностного аналога математической модели 
фильтрации используется консервативная абсолютно устойчивая неявная 


схема вычислений [2, 4, 5]: 


п+1 2 п п+12 п+12 п+12 п+12 
= в в-Ры — Реј ре Рају — Рр Реј Рл, 4: 
Е Г, 0.5 1+12,) Љ 1–1 2,ј Ах? 
| (7) 
п п п п 
ЈИ Руж Реј _ Рут Р,ј- 
на + 4: 1, ]+12 ла? 17-12 др? 
п+] п+1 2 п+1 2 п+1 2 п+і 2 п+1 2 
В, В РЕ; а Реј Ру __ Ру -РВ-ь) Т 

5 і, Ј 0.5 1+12,) Ах? 1-1 2,ј Дх? 

(8) 

п+1 п+1 п+1 п+1 
ера: < ЕТ СЕЈ = Рој — Реја ом 
г, /+12 Л? 1,3-1 2 2 7 
Функция Му ; моделирующая работу скважин, определяется 
следующим образом: 
~ ; | 1 К | | | 
Муза, Уа (9) | 
<. А ТА 2 
| Но 0А ур 


где 4; - дебит скважины в ячейке (положительный - нагнетательная, 


отрицательный - добывающая). | - 
Для определения насыщенности компонента в ячейке необходимо 
определить его расход через границы ячейки, который определяется явной 
разностной аппроксимацией формулы (4), а насыщенность компонента - 
аппроксимацией формулы (5). 
Для 


и граничных условий. 


решения уравнения (1) необходимо задание начальных · 


В соответствии со схемой движения жидкости между скважинами (см. 
рис. 1) на верхней и нижней границе рассматриваемой области может быть Е 


задано условие непротекания через. непроницаемые подошву и кровлю 


м лим 
а Дем 


> 


пласта. На левой и правой границе области также можно использовать 
такое условие, предполагая распространение потока только через 
рассматриваемую часть вертикального среза пласта. 

В этом случае граничное условие задается в виде: 

Ру а = Ру » (11) 
т.е. давление за границей области ру. равно давлению в граничной 
ячейке, что в соответствии с уравнением (4) дает нулевую скорость 
перетекания вещества через границу. 
~ В качестве начальных данных необходимо задать распределение 
давления и насыщенности компонентов в нулевой момент времени, их 
вязкость, режимы работы нагнетающей и добывающей скважин (дебит или 
забойное давление, возможную цикличность работы скважин), 
особенности вскрытия пласта в вертикальном разрезе, а также 
характеристики самого пласта (проницаемость и пористость). 

В приведенной модели используется модель неоднородного по 
проницаемости пласта с двумя видами пустотности (системой 
макротрещин и микротрещиноватостью блоков, что предполагает наличие 
двух различных значений пористости и проницаемости) (рис. 3). 

Для задания двух видов пустотности расчетная область покрывается 
сеткой, достаточно мелкой для того, чтобы задать различие свойств для 
матрицы и системы трещин. Структура матрицы моделируется 
объединением нескольких соседних ячеек с заданием для них одинаковых 
СВОЙСТВ. 

Ячейки вокруг матрицы играют роль системы трещин, для которых | 


характерны другие значения проницаемости и пористости. 


Модель пласта с двумя видами пустотности 
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а) б) 
а) Образец реальной структуры пласта (1 – трещина; 
2-матрица (блок); 3-каверны); 

6) Идеализированная модель (1-матрица; 2-трещина). 


Рис. 3 


Количество ячеек в матричной структуре (размеры матрицы 
коллектора) и различие свойств отдельных матриц между собой 
учитывается с помощью методов и соотношений теории вероятности на 
основании средних значений характеристик пласта и коэффициентов их 
вариации. 

На рис.4 представлены результаты тестового расчета вытеснения 
нефти из коллектора с одним видом пустотности (т.е. при отсутствии 


матриц). Стрелками показаны направления подачи воды и забора нефти. 


АИТ 
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3 ОПИСАНИЕ ПРОГРАММНОЙ РЕАЛИЗАЦИИ МОДЕЛИ 


Предложенная математическая модель реализована в виде 
исполняемого модула, предназначенного для работы на ІВМ-совместимых 
компьютерах в среде Мтдомљ 9х. 

Работу с программой можно разбить на несколько этапов: 

1. Задание размеров вертикального среза пласта и генерация сетки на 
расчетной области. 

2. Формирование внутренней структуры пласта (задание границ 
непроницаемого контура, блоков матриц и их размеры). Этот этап может 
выполняться либо визуальным заданием требуемых границ на расчетной 
сетке, либо автоматической генерации блоков матриц на основе 
указанных значений размеров и их вариаций. Пример расчетной сетки 


для некоторой гипотетической области после выполнения этого этапа 


показан на рис. 5 


Вид расчетной сетки 
для гипотетического пласта 


1-непроницаемые области; 2-блоки матриц; 3-трещинные области. 
Рис. 5 
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3. Задание особенностей вскрытия пласта нагнетательной и добывающей 
скважинами (на какую глубину, по всей толщине пласта и т.д.). 
Согласно расчетной сетке (см. рис. 5) крайние столбцы сетки 
представляют скважины, и программой предусмотрено указание любой 
ячейки из этих столбцов в качестве нагнетающей или добывающей 
(рис. 6). 


Особенности расположения 


нагнетающих и добывающих скважин 


ячейки; ||! -нагнетательная == - добывающая 
Рис. 6 


4. Задание режимов работы скважин - забойное РРСРР или дебит. 

5. Задание свойств пласта. Программой предусмотрены возможности 
задания двух видов проницаемости и пористости в соответствии с 
принятой моделью, анизотропии свойств для обоих видов структуры 
пласта, а также возможности вариации этих параметров в соответствии с 


различными законами распределения случайных величин. 


> 
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6. Задание условий окончания расчета - либо конечное время , либо. 
условие безводной добычи нефти. 

7. Задание периодичности вывода результатов расчета. С указанной 
периодичностью программа сохраняет значения давления 


и водонасыщенности для каждой ячейки расчетной схемы. 


Вт. Особенностью программной реализации данной модели является 


максимальное использование различных достаточно легко задавать 


и изменять исходные данные для расчета. 


результаты расчета полностью сохраняются в текстовом файле 
и в дальнейшем могут быть обработаны различными программами, 
например, Ехсе, Ѕш#ег. Отдельные графики могут быть построены 
непосредственно в процессе расчета, например, суммарный дебит 
скважины, пластовое давление в определенной ячейке. Это позволяет 
более оперативно оценивать результаты расчета и вносить коррективы 


в исходные данные. 
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